NACA RM No. 8107 


tl 


I v authority 


. . Restriction/Classification 
Cancelled 


Oupy l>Ju. f 

RM No. 8107 


* — 


CLASSIFICATION 
j NCLAbSfFiED 




c 


R E S EMjSfr W mS RAN D U M 

*CQ > \f***>' 

ihn t /o-C-s's.f^ 

RESPONSE OF A ROTATING PROPELLER 
TO AERODYNAMIC EXCITATION 
By 

Walter E. Arnold! 

Hamilton "Standard Propellers Division 
United Aircraft Corporation 
East Hartford, Connecticut 


Fo 


r Iterate* 


TCAsamED Docuvorr 

W B0T TO BS T.'JEN FRiM TIIIS 

*fT*=U^ Om Hattcml DsOaje jrf tb* United 
' — li* Act, 

i cr tte 


“^Restriction/ 

rralafcl: 


tofc^Classification 
“^'Cancelled 

'nvlliaa i 
G :t+ m r 

tK*r*Uu in up t-nuai aura cuue ai of knew* 
Vanity and dtacMttcn vtz =4 naewBlty eust b* 
wrfomact tbcracf. 


bjr Usr. 
Qftrtad 
d U 7 *l 
■Cprtkt* 
Fadsrtl 
later— t 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


WASHINGTON 

January 21, 1949 


Restriction/Classification 

*J&ai2sw. 


UNCLASSIFIED 


NACA EM No. 8107 



NATIONAL ADVTSQHT •COMHTi’l'EE FOE AERONAUTICS 


RESEARCH MEMORANDUM 


RESPONSE CF A ROTATING PROPELLER 
TO AERODYNAMIC EECIHAIION*' 

By- Walter E. Amoldl 

SUMMARY 


The flexural vibration of a rotating propeller "blade with clamped 
shank Is analyzed with the object of presenting, in matrix form, equations 
for the elastic "bending moments in forced vibration resulting from aero- 
dynamic forces applied at a fixed multiple of rotational speed. Matrix 
equations are also derived which define the critical speeds and mode 
shapes for any excitation order- and the relation between critical speed 
and blade angle. Eeference is given to standard works cn the numerical 
solution of matrix equations of the forms derived. 

The use of a segmented blade as an approximation to a continuous 
blade provides a simple means for obtaining the matrix solution from the 
integral equation of equilibrium, so that, in the numerical application 
of the method presented, the several matrix arrays of the basic physical 
characteristics of the propeller blade are of simple form, and their 
simplicity is preserved until, with the solution in sight, numerical 
manipulations well-known in matrix algebra yield the desired critical 
speeds and mode shapes frcaa which the vibration at any operating condition 
may be synthesized. 

A close correspondence between the familiar Stodola method and the 
matrix method is pointed out, ihdicating that any features of novelty 
are characteristic not of the analytical procedure but only of the 
abbreviation, condensation, and efficient organization of the numerical 
procedure made possible by the use of classical matrix theory. 


INTRODUCTION 


This report presents a theoretical analysis of the flexural vibra- 
tion of an aircraft propeller blade subjected to harmonic aerodynamic 
exciting forces at a fixed multiple of propeller rotational frequency. 

*This report Is a reproduction of Hamilton Standard Propellers* 
report No. HSP-613 , of December 11, 19^-7 , with some slight 
modifications to conform more nearly to NACA form. 
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Moat direct practical application is found in the calculation of response 
to rotational frequency excitation below the first mode resonant frequency, 
but the analysis also provides a means for studying critical speeds at 
any frequency order. However, since current interest is centered on the 
former application, coupling of blade bending and torsion is neglected as 
being of little importance when operating far below the first torsicmal 
natural frequency, and damping is likewise ignored. Initial offset and 
sweep are also neglected, although their effects cannot be dismissed as 
a generality. 

The propeller blade in this analysis is clamped at the shank, that 
is, constrained so that at some fixed location near the center of rotation 
only uniform rotational motion is possible. This condition of end fixity 
is chosen in accordance with the object of applying the results of the 
analysis to the particular case of vibration in response to air stream 
angularity, where a propeller of three or more bladeB vibrateB in an 
unsymmetrical mode which fulfills this condition. The "reactionlees modes" 
occurring in propellers of, four or more blades are also covered by this 
boundary condition. The Bystem is furthermore assumed to be linear, 
with small vibratory displacements, and simple bending theory is used, 
since the twist in convent i anal propeller blade designs is moderate. 

The differential equations of equilibrium are first derived and 
transformed into integral equations, and they are then examined in the 
form of a matrix equation for, a segmented blade, which permits the 
evaluation of vibratory response, critical speeds, and normal modes by 
simple classical methods which are particularly attractive because the 
solution in numerical form follows the symbolic form very closely. 


PHILOSOPHY OF ANALYSIS 


It must be recognized at the outset that the mathematical description 
of the vibration of a twisted, rotating, tapered beam, subjected to 
distributed vibratory loading, is essentially a complicated process, and 
that care must be taken to avoid dealing with expressions so cumbersome 
and involved in notation that the physical meanings of the various terms 
are completely hidden. The analysis here presented attempts to avoid 
such difficulties by using the concise abbreviations provided by simple 
matrix algebra, setting up matrix arrays of physical quantities wherever 
possible, and separating the operations of integration and differentiation 
from the physical quantities by the use of operational symbols in matrix 
form. It is thus possible to avoid writing large arrays of simultaneous 
equations, replete with multiple integrals, which might otherwise tend to 
disguise, through their complexity, the basic manipulations leading to a 
solution. By using concise matrix terminology, the problem is reduced to 
a classical form of equation which can easily be solved to yield Buch 
fundamental information as the critical speeds and normal modes of the 
vibrating system. As in simpler, ncnrotating systems, a knowledge of the 
normal modes and critical frequencies is sufficient to synthesize the 



NA.CA EM No. 8107 


3 


response of the system to an aerodynamic excitation, at any operating 
condition. 
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LIST OF SYMBOIS AND UNITS 

"blade station radius, inches 
"blade tip radius, inches 

"blade vibratory displacement, scalar, also matrix column, inches 

circular frequency of vibration, radians per second 

circular frequency of rotation, radians per second 

mass per unit length of "blade, also diagonal matrix, pound- 
seconds 2 per inch 2 

Young 1 s modulus, pounds per inch 2 

hlade section moment of inertia, scalar, also diagonal matrix, 
inches^ 

blade section angle from, plane of rotation, radians 
blade angle at reference station, radians 
twist matrix of sin 6 and cos Q 
section shear force, pounds 

section tending moment , scalar, also column matrix, pound- 
inches 

centrifugal tension, scalar, also d i agonal matrix, pounds 

applied air force per unit length, scalar, also column matrix, 
pounds per inch 

'integral operator, scalar, also matrix as defined later, inches 
derivative operator, scalar, also matrix as defined later, inched - -*- 

altered centrifugal tension matrix, -%T, pound-seconds 2 

O 

frequency matrix, seconds - ^ 
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P frequency order, cn/fl 

q)p frequency order matrix 

C curvature, scalar, also column, matrix, inches - 4 

n number of blade segments 

I unit matrix 

\ scalar frequency parameter in matrix equations, seconds - ® 

k modal column (moments), pound— inches 

k modal row (curvatures), inches - ■ 

u dynamic matrix, seconds - ® 

J orthogonal unit matrix, J 2 = —I 

Subscripts and matrix configurations are explained in the text and 
diagrams. Matrix notation follows the convent inns generally adopted in 
reference 1, where possible. 

EQUATIONS OF EQUILIBRIUM 


The forces and moments acting upon a differential blade element, 
in the plane of rotation, are shown by figure 1. The differential 
equations of equilibrium are found by separately summing horizontal 
forces, vertical forces, and moments from this diagram. Equilibrium of 
forces in the direction of the x-axis is expressed by 

dT + fi^px dx = 0 (l) 

and in the y-direction, by 

uFpy dx + dQ + F dx + ^^£} 2 px dx + = 0 (2) 

while the moments are summed by 


Q dx + dM = 0 


(3) 
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from, which may he written the three differential equations of 
e quillbr ium , 


dT p 

— + a PX = 

dx 


0 


00 


(a 2 + a£)py + H + 3P + = 0 


(5) 



dx 


( 6 ) 


Equation (if) may he integrated to obtain th6 centrifugal tension 
at any point along the blade 


r 


T = ft / px dx 
“x 


(7) 


Integrating equation (6) and substituting the expression for Q, 
derived by integrating equation ( 5 ) gives ■ 


M = (ft 2 + a£) / dx J~ pydx+ dx J I dx - / "S 


dx 


( 8 ) 


For vibration normal to the plane of rotation, the equilibrium 
diagram, of figure 1 would be altered only by the change In direction 
of the centrifugal force, ft^px dx, which would then be parallel to 
the x— axis. This eliminates the fourth term of equation ( 2 ), for 
vibration normal to the plane of rotation. Hie equation of moments 
is then 
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where it must be understood that the M, y, and F symbols represent 
different quantities from, those used in equation (8). 

Nov condense the notation by introducing the following operational 
symbols : 


Sf(x) 



f(x) dx 


Xf(x) = f*(x) 
dx 

Since equations (8) and (9) shall henceforth be used as simultaneous 
equations , introduce also subscripts p and r .to denote quantities 
referred, respectively, to vibrations normhl to the plane of rotation 
(parallel centrifugal field) and in the plane of rotation (radial 
centrifugal field). Equations (8) and (9), in reverse order, then 
became 


Mp = a£s2py p + S 2 F p - SXDy p 

(10) 

M,. = (n 2 + a^)S 2 py r + S 2 F r - 9TDy r 


Further to condense these equations, it is convenient to apply a 
matrix notation, in which 


Mpr 


{“4' 3pr = < 

fy-p] 

i i 

» F pr _ * 

<P = F 

Mr) 

[^rj 


frj L! 




n 


2 + o£ 



o . 


T 0 


s o' 


D 0 

P = 


n 


ll 

to 


, D = 



0 p_ 


0 T 

t • i 

r o 

01 


0 D 


The symbols, p , T, S, snfl D may be used either as matrices or as 
matrix elements without ambiguity, since the equations in which they 
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T rill "be used ■will make their meaning s evident. The matrix equation is 
thus 


Mpr - cpS^pyp r + S F^ STDy^pj, 


( 11 ) 


In order to solve this equation for tending moments in terms of applied 
air loads, it Is necessary first to eliminate the deflections, yp r , hy 
means of an additional moment— deflection relationship. 


FiaSHRAL EIGIDm OF A TWISTED BEAM 


Simple tending theory provides the scalar relation. 


M = EEC 


/ 

which applies to an untwisted team, where tending takes place atout one 
of the principal axes of inertia. This relation will te used in the 
twisted propeller tlade, tut it will te necessary to transform, the 
moment and curvature from, the twisted coordinates defined ty the principal 
axes to idle untwisted pr coord in ates. Denoting moments atout the 
minor axis .of inertia ty the sut script, f (for "flatwise" tending), and 
moments atout the ma^or axis ty e (for "edgewise" tending), the 
coordinate transformation may te written as follows, referring to the 
orientation of axes shown in figure 2. 

*%e = 0M pr 


where 


cos 0 —sin 0 


e = 


|^sin 0 cos 0j 

The curvature, expressed as a vector column matrix, may likewise te 
transformed in the same fashion — 


Cfe = ec pr 



8 


NACA EM No. 8107 


and these coordinate transformations may he introduced into the moment — 
curvature relation. 




9Mp r - EIg. e ©Cp r 

■whence 


C pr = ^(E 

But the curvature, in untwisted coordinates, is given by the second 
derivative of the deflection. 


c pr = ^pr 

Therefore 




The curvatures may be integrated to obtain deflections, using the 
operator. 


D X f (x ) = [ f (x) dx 
J 0 

which, although an integral operator, differs from the S already 
defined, due to the necessity for integrating curvature and slope from 
the origin outward. The desired relation between moment and deflection 
is then 


V 


r>r 





Pr 
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Substituting tliis into equation (ll) yields 


Mp r = qfi 2 pIT20-l(EI fe )-l0M pr + S 2 F pr - SUDD - 2 9~^- (e^ g ) ~^ r 


or 


Mp r = (<ps 2 pir2 _ anr- 1 )© 1 (xL fe )- 1 m pT + s 2 ^ ( 12 ) 


This is an integral equation in operational matrix foiUj •with, x as an 
independent variable, relating elastic bending moments to applied aero- 
dynamic forces. It can be solved by any of several methods, but this 
presentation will be confined to the discussion of a solution 
obtained by considering the propeller blade divided into a finite number 
of hinged segments. This treatment has the advantage of simplicity in 
directly reducing the matrix operators, S, D, and D -1 to convenient 
numerical form.. 


SEGMENTED BLADE 


Consider the application of equation (12) to a propeller blade 
divided into n equal segments, having its distributed mass divided 
proportionately among the hinge points and having par all lei springs 
across each hinge to represent the edgewise and flatwise flexibilities. 
Hie operators, S, D, and D - ^, can then be defined in series form., as 
follows : 


/ l n 

f (x) dx = > f (xj) ^ 

*i 


, d , , f (*i) - f O c i-i) 

“(a) * £ 


r xi i 

0 _1 f(xi) = / f (x) dx = ]>_ f (xj) Ax 

J o <3=1 
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The subscript nomenclature associated with these definitions is 
described by figure 3« Note that , because the moments in the segmented 
system are each one station removed from the forces which they equilibrate 
in accordance with equations (l), (2), and (3).» the subscripts for 
moments and inertias differ from the subscripts for masses, deflections, 
air forces, and centrifugal tensions at each station. 

There may be written n simultaneous equation's, corresponding to 
the n hinge points, in place of each of the two simultaneous equations (lQ). 
Instead of writing them separately, they will be combined into a single 
matrix equation by using the following conventions. 


MpOto) 

Mp 


y P (*i) 


faCxi) 


> , r p K, 1 -J 




> 


V^)J 


ypC*n) 




p! 0 - 0 


'0 

1 

0 


111-1 

0 

1 

w 

a 

0 

, T = 

0 t 2 - 0 

, S = Ax 

0 11-1 

— 


— _ _ 


0 0 1-1 



0 0 — T n 



_° 0 ~ p n_ 







0 0 0 -1 


1 0 0-00 


10 0-0 

-1 1 0-0 0 

3 

II 

110-0 

0-1 1-0 0 


111-0 

_ 0 0 0 1 1 


111-1 
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Similar definitions apply to the r— coordinates . Hote that K - ! is 
truly the matrix reciprocal of D. Equations ( 10 ) then are simultaneous 
matrix equations of rtr-order, and equation (12) will involve second— order 
matrices with n— order matrix elements, or, more concisely, equation (12) 
will "become a 2 n— order matrix equation. As an example, the twist matrix 
now "becomes 


cos 0 0 — 0 

—sin ©2. 0 — — 0 

0 cos Qg 0 — 0 

0 —sin ©2 — — 0 

0 0 cos S3 — 0 

— — — — — 

0 0 0 — cos 0^ 

0 0 — — —sin ©^ 

sin 6 j_ 0 — — 0 

cos 0-^ 0 — — 0 

0 sin Sg — — 0 

0 cos ©2 — — 0 

0 0 — — sin S rL 

0 0 — — cos © n 

— 

— 


and the stiffness matrix is 


XfCxo) 0 - - 

0 0 - - - 

1 O 

, * 

1 Ov 

t 
1 1 
1 1 

1 1 

0 0 - - - 

- - - 

— — — — — 

0 0 - - - 

Jr* 

V 

0 

1 
1 

1 

0 0 -r - - 

0 I e (xi) 


~ — — — -*-e (^h— l) 
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It is important to note that the reciprocals indicated in equation (12) 
are also of simple form, 0—1 being the transposed of 0 , and (Lf» e ) — 1 
being a diagonal matrix of the reciprocals of the individual inertia 
elements. Furthermore , I r~l is the transposed of S. 

Since , in propeller vibration studies, the bending moments referred 
to principal axes are directly related to measurable quantities, equa- 
tion (12) shall be rewritten in the fe coordinates, whence 

Mpg = 0 (cps 2 pir 2 - snr-lJe^tEIfe) -1 ^ + 9s20 _1 F fe ( 13 ) 


In solving this equation for .resonant conditions it will be necessary 
either to choose a fixed rotational speed and find the natural frequencies 
or to choose a propeller vibration order (multiple of rotational speed) 
and find the critical speeds. Adopting the latter procedure and 
letting P = write 



T = ft^r 


Also abbreviating the air moment column, 

= 0S 2 0- 1 F fe 

gives 


Mfe = fl 2 0(q)pS 2 pF- e - ST]T-l)8-l(EI fe )“ 1 Mf e + Mfe^ 
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Solving for Mf e 

Mf e = (1 - fl 2 e(cp p S 2 p3)- e - STEr 1 >- 1 (EX fQ )- 1 ) _1 M fe (a) 

lotting X = 
and 


u = 8(cppS 2 pD" 2 - SrIT- 1 )9- 1 (EE fe )- 1 


CUt) 


Kfe = X(XI - u)-\( a) 


(15) 


Numerical quantities can. be substituted into equation (15 ), and 
the bending moments caused by applied air moments at a predetermined 
rotational speed frequency order may be directly calculated. 
Another form of solution, obtained by the use of Sylvester *s Theorem, 
will often be more useful. This theorem, -which provides a series 
expansion of a matrix polynomial, applied to equation ( 15 ) yields 


Mfe = 


y 5 - A M (r>* 

ja ] X “ Xj» 


( 16 ) 


where M'- J is a normal mode moment component of 
terms of the modal row and colum n associated with 
from the relations. 


%e (a) , 
a latent 


obtained in 
root, Xp, 


(l r I — u)k r = 0 

(IT) 

— u) = 0 

( 18 ) 


*The subscript r is used henceforth to denote any of the 2n modes of 
the segmented system. This preserves similarity with matrix notation 
in reference 1 and should cause no ambiguity, s inc e the pr coordinates 
do not appear again in the report. 
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= (KjJCp) 1 k r K r %e^ a ^ 


(19) 


Equations (17) and (18) are known, alternatively as the "charac- 
teristic equations" of a matrix, u, and are typical of problems in 
■which there is an independent parameter, such as vibration frequency, 
as veil as -a group of coordinate dependent variables. The "latent roots," 
or values of the parameter for which there is a solution, and the 
coordinate values associated with each of these roots are frequently the 
solution desired, hence many classical methods are available for 
numerical solutions. In the case of the vibrating propeller blade 
represented by the characteristic equation ( 17 ), the latent roots of 
the matrix, u, are critical speeds (strictly, reciprocals of critical 
speeds squared) and the modal columns, kj., are the bending-moment 

distributions corresponding to the natural modes for the critical speeds. 
The modal rows, K r , defined by equation (l 8 ), are the corresponding 
curvature distributions, as will be shown. 

Numerical methods for determining the modal rows, columns, and 
roots are described in detail in reference 1. In the case at hand, the 
labor involved can be reduced by deriving a relation between the modal 
row and column, first forming the transposed of equation ( 17 ) . 


^•(Xrl -u») = 0 


Equation (l4) stated. 


u = e(cp p s 2 pir e - STD - 1 ) 9 -1 (life )- 1 

Since 


9* = 0" 1 , S* = B - 1, Tip* = <p p , p* = p, t» = t 

and 
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it follows that 


u‘ 


= (Ei fe ) _1 e(0 p s 2 pir^ - stirrer 1 


= (ECfe) ^fee) 


whence 


V(jx* - (EXfJ-^CEEfe)) = 0 


and 

k r * (EX fe )- 1 (X r I - u)(EX fe ) = 0 
or 


h r * (Elf e ) -1 Or I - u) = O 


Thus, the modal row is related to the modal column, by the equation, 

K r = kp^EXfe)- 1 (20) 

which shows that the modal row is a curvature distribution, since the 
modal column, from which it is here derived, satisfies the equation for 
the elastic moment. 


EFFECT OF BLADE MGLE OR CRITICAL SPEEDS 


For small changes in blade angle, the changes in critical speeds 
can be found without recalculating u for a new blade angle. A partial 
derivative of the characteristic equation ( 17 ) may be formed as follows: 


(>*I “ u)k r = 0 





o 
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Preimiltiplying "by K r , the second term vanishes, leaving 


■whence 


dp v ^ r dp 


( 21 ) 


New examine u and. find its derivative. 


u = ©(qa^pir 2 - STD -1 )er 1 (Ei fe )- 1 


H= | (V^ - ST D -l)9-l( E%e )-l 


+ ©((PpS 2 pD-2 - STEr-l)^|-i (Elfe ) -1 


( 22 ) 


Since 


© = 


cos 0 
sin 0 


—sin 9 
cos 0 


and d0 = dp 



cos 9 
—sin 9 


sin 0 
cos 0 


d© 

dp 


— sin 0 

—cos 0 

d© 1 

—sin 0 

cos 0 

cos 0 

-sin 0 

' SP 

-cos 0 

-sin 0 
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Let 


J = 




then 




5e 1 Q -i T 

1 =0 J 


Substituting these derivatives into equation (22), 


~ + uCEXfJjC^e) 1 


Premultiply by tip and postmultiply by k^, so that 


K r |^r = ~ K r Juk r + B r u C S ^e) J ( EI fe) "Sr 
and by meting use of the cbaracteristic equation again, u is eliminated. 


K r = -ttpJXjJs:,. + (Elf e ) J^KLf Q ) 1 k r 


The second team can be simplified by noting, from, equation (20), 

that 

V = (EIfe)~ lfc r 

and 


^r* = KpCEXfe) 
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■whence 


«r *r = + V J *rO 

But, since this Is a scalar equation, the elements in any term, may be 
transposed, so that 


and 


kj* * Jk p 



kj. = -SXyfCjJkj. 


whence, equation (21 ) becomes 


1,e r Jl£ r < 2 3) 

As a last refinement, note that 

1 ^ fi r 

dp dp dp 

r 


dn r 2 

dp 


1 ^ 

d2 





n r dp 


= (K-rkr) 1 K r Jk r 


(2b) 


This relationship makes it possible to make use of a solution, 
obtained for any arbitrary blade angle, to find the nature of the varia- 
tion in critical speed with variations in blade angle, at least for a 
small range in angle. It is interesting to put this equation into an 
integral form, for the distributed system, by noting the meanings of the 
operations indicated by equation (2h). The scalar, Kjk r , is a sum of 
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the products of each pair of curvatures and moments and hence represents 
twice the strain energy in a natural mode, so that, in terms of scalar 
moment elements, Mf and Mg, there may he written 




dx + 



dx 


Furthermore, the operator, J, has the effect of interchanging flatwise 
and edgewise (f and e subscript) quantities, so that 


KpJtp 5 




^e 


dx 


and the changes in critical speed may therefore he expressed as a function 
of a small change in blade angle hy the following scalar equation of 
definite integrals evaluated at the critical speed: 


«r 



a-P 


C25) 


This useful relationship might have been obtained by other methods, 
but the presentation here given has seemed the most straightforward to 
the writer. 


NUMERICAL COMPUTATIONS 


The foregoing analysis has provided several equations in matrix form 
which may be employed directly in numerical work. After forming the 
basic arrays of physical quantities and combining them to obtain the 

matrix, u, and the column, equation ( 15 ) may be solved for the 

vibratory bending moments at the n blade stations by methods such as 
proposed by Aitkan (reference 2) in order to find the nonre sonant response 
of the propeller to the applied air forces at a given frequency order. 

If the critical speeds are desired, the characteristic equation (17), may 
he solved hy repeated premultiplication of' an arbitrary col umn , as 
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described in reference 1 , and "by this means the bending-moment distribu- 
tion of the normal mode is also obtained. The normal mode curvature 
distribution may be found by similar poatmultiplication of an arbitrary 
row, according to the characteristic equation ( 18 ), or directly from the 
moment column, by equation (20). The critical speeds, normal, mode 
columns and rows, and the applied air moments may be gathered together in 
equation (l 6 ) in order to express the resultant vibratory bending moments 
as functions of rotational speed. 

For numerical purposes, it is often conven i ent to expand equation ( 16 ) 
in a form which will converge more rapidly, thus requiring the calculation 
of fewer mode shapes. It is easily shown that 


%e = % e ta) + 2 


M^ r) 


r=! x ~ hr 


which follows from observing that 


(26) 


lim 
ft 2 — >0 

(x, 


2 n 

“t» - 2 . 


r=l 




Since the matrix, u, includes an arbitrary blade angle, the 
numerical solution in terms of normal modes and frequencies yields 
complete information at a given frequency order only for a single blade 
angle,* however, for smal 1 variations in blade angle, equation (24) may 
be used to determine the effect on critical speeds without going through 
a completely new solution. 


In computing critical speeds for a fixed vibration order, it . should 
be noted that same of the 2n latent roots, Xp, may be negative. In the 
particular case of first-order vibrations, not more than one of these 
roots will, in general, be positive. The occurrence of negative roots 
implies imaginary critical speeds, which are difficult to picture as such, 
but may appear more logical if treated as the positive critical speeds 
which would exist were the centrifugal forces and inertia forces to be 
reversed. This condition is probably more easily understandable for the 
special case of P = 0 , when the negative roots correspond to critical 
speeds for several modes of buckling under reversed centrifugal forces. 
Figure 4 presents a sketch of the natural frequency spectrum of a propeller 
as a function of rotational speed, using frequency squared along each 
axis, so that the negative roots appear as intersections in the third 
quadrant . 
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STODOLA* S METHOD 


la tlie calculation of natural frequencies in problems Involving tlie 
Titration of teams, the Stodola process is often employed . This is 
normally a method, by 'which, a deflection curve is assumed and used to 
calculate the inertia loading, which, is then successively integrated 
along the length, of the team in order to ottain the shear, moment, slope, 
and a new deflection shape. The new deflection curve is then carried 
through the same process again, and repeated iterations final l y converge 
on the true mode shape, with an increase in amplitude through each itera- 
tion in proportion to the lowest natural frequency. After solving for 
the lowest mode and frequency, higher modes can also be found, hut since 
there is a tendency for the fundamental mode to became more prominent 
through each iteration, it is necessary to employ the orthogonality 
relations in order to eliminate the unwanted mode or modes in each 
higher mode solution. 

The Stodola process could sIbo be used with an assumed moment distribu- 
tion, carrying through each series of integrations to obtain an improved 
moment mode shape, and it is interest ing to note that the solution of the 
matrix equation (17 ) for its latent roots and modal columns is very closely 
related to this type of Stodola process. Premultiplication of an 
arbitrary moment column, by u, upon examination of the definition of u, 
equation (l4), includes in a lumped form the several steps of successively 
integrating the curvature, slope, loading, and shear, including centrif- 
ugal tension, effects, to obtain an improved moment distribution, and the 
comparison of successive mode shapes, after convergence has been attained, 
yields the lowest natural frequency in the form of the dominant latent 
root. The process of modifying the u-matrix in order to eliminate the 
dominant mode and permit solution for the sub dominant root, when care- 
fully examined is found to provide a new matrix which represents the same 
series of operations during each iteration, with the addition of an 
"orthogonal! zing" step Included in each iteration. TJse of the Stodola 
process in this manner has long since been successfully accomplished at 
least in the somewhat simplified case of a nonrotating untwisted beam 
(reference 3)- The rotating, twisted propeller blade requires numerical 
calculations of larger volume and hence greater difficulty, but no new 
principles are involved. 

The matrix method of calculation, therefore, is not basically novel, 
but represents a technique for organizing in an efficient manner the 
multitude of operations involved in a large-scale numerical calculation. 

One of its strangest attractions is that it permits a proper perspective 
to be maintained not only in the analytical p ha ses of the problem but 
also throughout the numerical work. 
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Figure 2.- Orientation of twisted (fe) coordinates in relation to unt wisted 

(pr) coordinates. 
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Figure 4.- 




